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We present an operational definition of the Wigner function. Our method relies on the Fresnel 
transform of measured Rabi oscillations and applies to motional states of trapped atoms as well as 
to field states in cavities. We illustrate this technique using data from recent experiments in ion 
traps [D. M. Meekhof et al., Phys. Rev. Lett. 76, 1796 (1996)] and in cavity QED [B. Varcoe et 
al., Nature 403, 743 (2000)]. The values of the Wigner functions of the underlying states at the 
origin of phase space are Wjo)(0) = +1.75 for the vibrational ground state and W|i)(0) = —1.4 for 
the one-photon number state. We generalize this method to wave packets in arbitrary potentials. 



Currently the Wigner function [lj enjoys a renaissance 
in many branches of physics, ranging from quantum op- 
tics 

2], 

nuclear Q and solid state physics to quantum 
chaos |4J. This renewed interest has triggered a search 
for operational definitions [f| of the Wigner function, that 
is, definitions which are ba7ed on experimental setups @. 
In the present paper we propose such a definition based 
on the resonant interaction of an atom with a single mode 
of the electromagnetic field. In contrast to earlier work it 
is the atomic dynamics that performs the major part of 
the reconstruction. Our definition is well-suited for the 
Jaynes-Cummings model p| but allows a generalization 
to other quantum systems. 

Many methods to reconstruct the Wigner function 
of a cavity field or the motional state of a harmonic oscil- 
lator have been proposed jg]. In particular, the Wigner 
function of a cavity field can be expressed in terms of 
the measured atomic inversion [jj. This operational 
scheme 01 lives off the dispersive interaction between 
the atom and the field. Since this scheme requires long 
interaction times, an operational definition of the Wigner 
function based on a resonant interaction is desirable. The 
method of nonlinear homodyning 11] and quantum state 
endoscopy 0] fulfill this need, but require rather com- 
plicated reconstruction schemes. 

In contrast, the Fresnel representation proposed in the 
present paper is rather elementary. It expresses the value 
of the Wigner function at a phase space point a as a 
weighted time integral of the measured atomic dynamics 
caused by the state displaced by a. The weight function 
is the Fresnel phase factor. Since this method can be 
applied to a cavity field as well as a trapped ion |13| we 
use a harmonic oscillator as a model system. However, 
we show that this approach can easily be generalized to 
non-harmonic oscillators. 

Our definition relies on controlled displacements pH 
IhH of the quantum state of interest and the observation 
of Rabi oscillations of a two-level atom interacting reso- 
nantly with this field. We record the probability 
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of finding the atom in the ground state \g) as a function 
of dimensionless interaction time r and complex-valued 
displacement a. Here P n {ct) = (n\D(a)p D^{a)\n) de- 
notes the occupation statistics of the state p displaced 
by the displacement operator D(a). 

The Wigner function VF(a) of the original state p is 
determined ^(| by the alternating sum 
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of the probabilities P n . Hence, we can find the Wigner 
function from the atomic dynamics when we measure the 
probabilities P g for different interaction times and solve 
Eq. for the occupation probability P n {a). 

However, there is no need to evaluate the probabilities 
P n explicitly ^tJ- We can obtain the Wigner function 
directly from the atomic dynamics P g without ever cal- 
culating P n , making use of the integral relation 
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When we multiply P g — 1/2 from Eq. by 
2/(7r-\A)exp [ir 2 /it] and integrate over r, Eq. © yields 



Pg(r;-a)-- 



1 OO 

= -£(_l)»P n (-a). 



We recall the connection, Eq. J5J), between W(a) 
and P n (a) and find the Wigner function VF(a) = 
lim W(a; r rn ) in terms of the truncated Fresnel trans- 
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of the Rabi oscillations. 

Hence, the Wigner function at the phase space point a 
is determined by a weighted time integral of the atomic 
dynamics due to the initial state displaced by —a. The 
weight function is the Fresnel phase factor exp(ir 2 /7r). 
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FIG. 1: Fresnel representation of Wigner function illustrated 
by the first excited energy eigenstate of a harmonic oscillator. 
The atomic dynamics P g (r;a), Eq. JQ, as a function of in- 
teraction time t and real-valued displacement a (top) serves 
as input to the truncated Fresnel representation, Eq. @. 
For a finite measurement time r m the so-obtained function 
W(a; Tm) is complex-valued. The real part of W(a; r m ) (bot- 
tom) approaches the correct Wigner function, W^, already 
for moderate measurement times. 




FIG. 2: Asymptotic approach of the truncated Fresnel repre- 
sentation for the first excited state. Real part (black line) and 
imaginary part (gray line) of W(a = 0;r m ) oscillate around 
their asymptotic values —2 and 0, respectively, with slowly de- 
creasing amplitudes. The asymptotic values emerge already 
after a measurement time of the order of 2n. 



It is instructive to compare the Fresnel re pre sentation, 
Eq. i to the Fourier method used in Ref. [l3| . The lat- 
ter relies on the analysis of the Fourier transform of the 



Rabi oscillations P g and thus needs the full functional 
dependence on frequency. In contrast, the Fresnel repre- 
sentation gives the Wigner function as a single integral of 
P g without the need to analyze an intermediate function. 

The Fresnel representation seems to suffer from three 
disadvantages: (i) It contains the complete time evolu- 
tion of P g , that is, from r = to r = oo. However, 
any experiment can only record the dynamics for a finite 
measurement time r m . (ii) It can take on complex val- 
ues. However, the Wigner function is always real. (Hi) 
It relies on a continuous time evolution. However, the 
experiment can only provide a discrete sampling. 

We address each of these problems separately and 
first demonstrate that the phase factor exp(ir 2 /7r) makes 
the integral insensitive to the long time behavior of P g . 
For this purpose we show in Fig. the Fresnel recon- 
struction for the first excited energy eigenstate of a 
harmonic oscillator. The corresponding Wigner func- 
tion W\x){a) = -2 (l- 4|a| 2 ) exp(-2|a| 2 ) is rotation- 
ally symmetric. Therefore, it is sufficient to depict W\\\ 
along the real axis. 

The top of Fig. n displays the atomic dynamics P g as 
a function of interaction time r and real- valued displace- 
ment a. Sinusoidal Rabi oscillations appear along a = 0. 
Moreover, for appropriately large displacements we find 
collapses and revivals. In the bottom part we present the 
real part of W(a;r m ). The time axis corresponds to the 
upper limit r m of the truncated Fresnel representation, 
that is the measurement time. For very short times the 
function has no similarity with the correct Wigner func- 
tion, W\i\. However, for longer measurement times the 
curves approach the exact Wigner function. 

In order to study the asymptotic behavior in more de- 
tail we show in Fig. [21 a cut of the truncated Fresnel rep- 
resentation along a = 0. Here we depict both, real and 
imaginary part ofW(a;T m ). We note that after r m ~ 2-7T 
the real part, indicated by the black line, begins to os- 
cillate around the correct value 1 / F|i)(a = 0) = —2 with 
decreasing amplitude. In a similar way, the imaginary 
part, denoted by the gray line, approaches its asymptotic 
value zero. Unfortunately, this convergence is slow. This 
feature is a consequence of the familiar Cornu spiral 0, 
also used by R. P. Feynman and J. A. Wheeler to ad d up 
scattered waves — a precursor of the path integral |19j . 
Despite the slow asymptotics the knowledge of the Rabi 
oscillations up to a measurement time of a least r m = 2tt 
is sufficient to obtain the asymptotic Wigner function. 

The upper part of Fig. [3| shows measured Rabi oscilla- 
tions of an ion due to its center-of-mass motion [2fl l2l| . 
Due to the high sampling frequency we can directly inter- 
polate the measured data to evaluate the Wigner func- 
tion. In the bottom part we display the approaches of the 
real part (black line) and imaginary part (gray line) of 
the truncated Fresnel integral towards their asymptotic 
values (dashed lines). For the ion in the motional ground 
state with Wigner function W| )(a) = 2 exp(— 2|a| 2 ) 
we find the asymptotic value VF|o)(0) = +1.75 + i 0.05 
whereas the ideal value is VF|o)(0) = 2. The presence 
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FIG. 4: Dependence of the Wigner function W^(a = 0) at 
the origin of phase space on the cutoff TV of the overdeter- 
mined system, Eq. @. The values scatter around VT|i)(0) = 
— 1.4. The Rabi oscillations |23j| (inset) result from a cavity 
field prepared in a one-photon number state. 



FIG. 3: Wigner function at the origin of phase space obtained 
from experimental Rabi oscillations (top) of a stored ion |2(i| 
using the truncated Fresnel representation, Eq. @: Approach 
(bottom) of real (black line) and imaginary parts (gray line) of 
W(a — 0;r m ) towards their asymptotic values (dashed line) 
yielding Wi o \(0) = 1.75-M0.05. The ion was prepared initially 
in the atomic and vibrational ground state. 



of an imaginary part in the estimated Wigner function 
does not represent a fundamental drawback, as long as 
it comes from a complex integral of realistic data that 
ideally should produce a real value. 

It is quite remarkable that the method works so 
well considering the fact that in the particular experi- 
ment |20| the atomic dynamics is g overned by the non- 
linear Jaynes-Cummings model |22|. In this case the re- 
lation between the atomic population and the photon 
statistics is not of the form of Eq. QJ. However, for 
low excitations the deviations from the familiar Jaynes- 
Cummings model are small as shown in Fig. 1(b) of 
Ref. [2(|. Therefore, superpositions of low excitations 
can be successfully reconstructed. 

The inset of Fig. |U shows Rabi oscillations of an atom 
due to a cavity field 23] . Since in this case only a few data 
points are available we have to employ a discrete version 
of the Fresnel representation. The data [23 consist of a 
finite set of M interaction times Tj and measured proba- 
bilities P g (rj\ —a). We now choose kernel coefficients fj 



such that 



^cos(2V^TTr J )/ 3 = (-ir +1 , 



(5) 



multiply Eq. Q evaluated at discrete times Tj by fj , and 
sum over j. With the help of Eq. JSJ we then arrive at 
the discrete representation 
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of the Wigner function. 

Since the photon number n runs from zero to infinity, 
Eq. describes a system of infinitely many equations 
for M unknown coefficients fj. When we consider the 
truncated system of N > M equations we find the so- 
lutions fj of this overdetermined system by minimiz- 
ing the sum of quadratic deviations. Substitution of the 
so-calculated coefficients /j into Eq. (JSJ results in the 
Wigner function W^ N '(a). 

Figure 0| displays the dependence of W^ N \a = 0) on 
the cutoff N. The values scatter around W / |i)(0) = —1.4 
indicated by the dashed line. We recall that the ideal 
value of the Wigner function corresponding to a one- 
photon number state is W w (0) = -2. For N = M = 10 
the method is extremely sensitive to uncertainties in the 
measured data resulting in a large deviation from -1.4. 

The error bars SPj translate into a standard deviation 
1 1/2 
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of the Wigner function. Accord- 



ing to Eq. © the coefficients fj are solely determined 
by the set of Tj and can take almost any values. Obvi- 
ously large values give rise to a large error in the Wigner 
function, unless the corresponding uncertainties SPj are 
small and compensate this effect. 

For the experiment of Ref. [23) performed for a differ- 
ent purpose and therefore not taking advantage of this 
feature, we find 5W ss 4.3. However, a future experiment 
can either choose the interaction times Tj such as to min- 
imize the fj or increase the accuracy of the measurement 
of P g at times where fj is large. 

The Fresnel representation, Eq. of the Wigner 
function as well as the discrete version, Eq. ©, rely 
on the resonant Jaynes-Cummings model. However, the 
technique to obtain the Wigner function directly by an 
appropriate integral transform of the time dependence 
of a measurable quantity is much more general. Indeed, 
our method even allows us to obtain the Wigner func- 
tion of a wave packet |^(0)) = Yin VVil^n) consisting 
of a superposition of energy eigenstates \<j> n ) with en- 
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ergy E n = hw n of a potential V = V(r). Here, we 
do not restrict ourselves to a harmonic oscillator po- 
tential. Our only assumption is that V is symmetric 
with respect to the origin. In this case, the represen- 
tation Eq. J5J of the Wigner function still holds true. 
Moreover, we recall that the autocorrelation function 
C(t) = (ip(0)\ip(t)) = J2 n I^Vil 2 exp(-icj„*) is measured 
routinely in wave packet experiments, e.g. j24[. 

When we apply the displacement D to the initial 
state p = \ip(0))(ip(0)\, the autocorrelation function 
C(t; a) = P n (ct) e~xp(—iui n t) is similar to Eq. with 
P n (a) = (<f> n \D(a)pD^ (a)\(p n ) . We multiply both sides 
of the equation for C by 
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: cos[n(w)7r] 



(7) 



and integrate over t. Here n(u>) is an arbitrary and con- 
tinous function such that n(uj n ) = n. 

When we make use of the symmetry relations f{—t) = 
/*(*) and C(—t) = C*(t), we find the Wigner representa- 
tion 



W(a) =4 Re 



dtf(t)C(t;-a) 



(8) 



of the initial wave packet in terms of the measured auto- 



correlation function C and the integral kernel /. 

Equation (JHJ) is the generalization of the Fresnel repre- 
sentation, Eq. (0J, to a quantum system with arbitrary 
discrete spectrum. For wave packets in a harmonic os- 
cillator or in a box with energy spectra cu n = 27m/T c ; 
or uj n = 2irn 2 /T r , respectively, the integral kernels from 
Eq. (7J are delta functions at half of the classical period 
T c i, or half of the revival time T r . For these potentials, 
the time evolution itself creates the function (—1)™, that 
is the parity operator. For a spectrum uj n = y/n + 1 2f2, 
Eq. J7J predicts a quadratic Fresnel-like phase. 

Complications arise for spectra with degeneracies, such 
as the hydrogen atom. Nevertheless, we can recon- 
struct Rydberg wave packets measuring the autocorrela- 
tion function at half of the revival time. These and other 
relevant examples will be discussed elsewhere. They sug- 
gest a wide application of integral (Fresnel) transforms 
in the search for a practical reconstruction of the Wigner 
function. 
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